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ABSTRACT 


Calculation  of  for  satellites  using  accommodation  coefficients  is  reviewed.  A 
idtenomenological  model  for  acconunodation  coefficients  due  to  Hurlbut,  Sherman,  and  NociUa  is  used 
to  obtain  values  for  the  accommodation  coefficients  for  average  satellite  materials,  thennosi^bere 
oonstitueots  and  tsmperatuies,  and  satellite  velocities  using  a  number  of  laboratory  measurements.  There 
is  a  significant  difTereoce  betweeni  these  results  and  the  traditional  method  of  calculating  These 
differences  contribute  as  much  as  20%  error  in  use  of  tbermospbere  models  for  calculation  of  satellite 
drag. 
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1.  INTRODUCTION 


The  drag  force  per  vinit  mass  (AO  on  an  element  of  area  (A)  of  a  satellite  is  given  by 


where  /)  is  the  atmospheric  density  and  V,  is  the  speed  of  the  satellite  with  respect  to  the  atmosphere. 
Analysis  of  satellite  drag  involves  the  product  of  three  quantities:  the  projected  area  to  mass  ratio  {A/M), 
the  ballistic  coefficient  (Q,  and  the  atmospheric  density  (/>).  None  of  these  is  known  without  error. 
General  analyses  {1-3]  suggest  that  the  error  in  calculation  of  drag  ranges  from  15%  at  3(X)-  to  600-km 
altitudes  to  30%  at  800-  to  1200-km  altitudes.  This  is  true,  even  for  spherical  satellites  where  the 
I»ojccted  area  is  known  without  error.  Due  to  the  complexity  of  thermosphere  models,  it  has  been 
generally  assumed  that  the  error  is  due  to  inadequacies  of  these  models.  There  is  abundant  evidence  that 
thermosphere  models  have  the.«»  errors  (2,3].  The  purpose  of  this  discussion  is  to  evaluate  the  possible 
error  duh  to  the  uncextalnty  in  the  knowledge  of 

There  are  some  limits  on  the  focus  of  this  report.  Cook  (4]  provides  an  excellent  review  of  the 
issues,  which  are  summarized  here.  The  mean  free  path  in  the  thermo^jhere  is  over  200  m  at  an 
alUhide  of  200  km  and  Inoeases  to  more  than  600  m  at  250  km.  If  the  satellite  linear  dimension  is  L, 
then  the  Knudsen  number  XJL»\  and  free  molecular  flow  occurs.  Therefcae,  consideration  is  limited 
to  calculation  of  for  free  mol(»»lar  flow. 

In  low  earth  orbit  (LEO),  the  kin^e  energy  of  molecules  relative  to  a  satellite  surface  is  about 
0.323  eV/aanu  or  about  5  cV  for  atomic  oxygen.  LEO  velocity  Is  aj^ximsaely  7.5  km/scc  for  circular 
orbits,  and  high  eccentricity  orbits  can  exceed  this.  Gwsynduonous  orbit  (GEO)  velocity  is 
approximately  3.0  km/sec.  From  Bguie  1  the  intere^  in  molecule  kinetic  energy  is  determined  to  range 
from  0.10  to  12  eV. 


0.1  1  10  100 
UOLECULE IQHET1C  imERQY  (oV) 

Figure  1.  Veiociiy  vs.  kineiic  energy. 
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Toe  thermosphere  constituents  of  concern  here  are  H,  He,  N,  O,  N^,  O2,  and  Ar.  At  200  km, 
and  0  are  the  major  constituents  and  have  nearly  the  same  number  density.  Above  600-km  altitude  H  and 
He  number  densities  are  conq)arable  with  O,  aud  Nj  has  decreased.  At  higher  altitudes  H  becomes  the 
dominant  constituent  These  results  are  shown  in  Rgure  2,  calculated  from  Hedin’s  MSIS84  thermosphae 
model  [4], 


Finally,  the  numerical  value  of  Q  is  differeifr  for  each  constituent.  Therefore,  the  drag  force 
should  be  written 
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Figure  2.  Moiecule  number  density. 


where  p,  is  the  density  for  the  constituent  is  die  (xmrespooding  ballistic  coeffidcnt 


2.  ACCOMMODATION  COEFFICIEP^rS 


The  general  theory  of  the  interaction  of  a  surface  and  a  gas  is  not  yet  at  the  state  where  a  physics- 
based  model  can  predict  die  exchange  of  energy  and  momentum.  Therefore,  one  must  make  use  of 
macroscopic  averages.  The  formulation  defines  three  accommodation  coefficients  defined  by  the  thermal 
accommodation  coefficient 


a=- 


ErE^ 


the  tangential  momentum  accommodation  coefficient 


(r=- 


and  die  normal  momentum  accommodation  oiefQcient 


PrPw 

£(  and  £,  am  the  energy  fluxes  Inchkni  on  and  reemiued  from  a  surface  clcmem  per  unit  time.  Ihc 
quantity  £.  is  the  energy  flux  that  would  be  carried  away  if  all  Uie  acconunodaied  molecules  were 
reemitted  with  a  Maxwellian  tfistiibuUon  with  temperature  T,  (the  surface  temper^ure).  So,  a  is  a 
measure  of  the  aowunt  that  molecules  have  their  energy  "acconiraodaticd"  to  whm  it  would  be  if  all  those 
reemmitted  had  an  energy  Rtr  complete  acccnimcdation  ce?:!,  while  whli  no  accommodation  op^O. 
Complete  acconmiodation  (o^l)  is  often  referred  to  as  diffuse  splattering,  whereas  the  latter  (oMl)  is  called 
scaUaing. 

The  treatnumt  of  momentum  transfer  origtoally  due  to  Maxwell  postulated  Uiat  a  fraction  (t-o) 
of  tlie  incident  molecules  bad  specular  scattering,  and  the  remaining  incident  fraction  (o)  was  scattered 
diflusely.  This  formulation  proved  to  be  inadequate,  which  led  to  intnxhicing  two  coefftcienfs  as  defined 
above  in  analogy  to  (he  expression  for  ol  in  this  case  t  and  p  arc  the  tangential  and  notmal  momentum 
components.  Now,  p«  and  are  the  norma!  and  tangential  momentum  components  romult^  with  a 
Maxwellian  distribution  Therefore.  x.=0.  For  the  Umiting  cases  of  entirely  specular  reflection— with 
no  energy  exchange— and  for  completely  diffuse  sc^teting.  . 

In  gcttuaJ,  these  throe  i»rameters  are  considered  independent  and  will  depend  on  prt^ties  of 
the  gas  and  surface  material  states.  In  addition  thete  is  considerable  evidence  that  these  parameters  vary 
with  iocidcni  angle.  Tbetefore,  (he  accommodaL..a  cocfTscienis  arc  considered  to  be  of  the  form 


o(0)“ 


m-m 
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and 


a(0)= 


riieyrxe) 

rX») 


a'iey 


Pi0)'Pw 


Tliese  accommodation  coefficients  are  quite  arbitrary,  and  a  number  of  authors  have  introduced  variants 
of  the  momentum  coefficients  that  dismay  their  data  to  advantage.  The  first  variant  [6]  defines  the 
"relative  niOTmal  momentum  fiansfer  RKT  as 


r  PrPr 
“r 


Pi 


Tim  second  variant  is  to  view  the  normal  momentum  coefficient  as  a  vector  quantity  and  define 


n. 


W 

Pl*Pm 


The  third  {E,9J  uses  both  these  ideas  to  define 


In  the  Itdlowing,  this  last  definition  is  tised  wHh  the  label 


In  all  (sses.  these  cocfildcnts  aie  assumed  to  depend  on  the  incident  angle  0. 

Schmnbag  (10,1 1}  introduced  (hscussion  of  the  senttoing  charactefistics  in  addihon  to  the  amount 
of  accommodation.  He  oon^dered  the  angle  of  the  (mean)  reemitted  beam  and  ttm  shaqjc  of  the  beam 
*V;  sec  Figure  3.  Huiibut  (8,12)  reviews  more  modem  data  on  these  pn^jciiies.  The  issues  remain 
Sdrambetg  discussed  a  general  distribution  fimetion  of  scattered  moieties  about  the  mean  reemitted 
beam  Because  genera!  models  were  not  avail^le  at  that  time,  be  was  only  able  to  satisfactorily  treat  two 
limiting  cases.  He  did  (his  by  intmdudng  a  malhemaii^  relation  between  the  &  and  9, 
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cos6/={cos9y 


231741-3 


Z 

Figure  3.  Flotation  for  surface  interaction  model 


Rrst  is  the  traditional  case  of  specular  reflection  for  which  This  is  obtained  when  y=1>  Here,  the 
particle  reflects  off  the  surface.  Second  is  the  case  of  diffuse  reemission  In  this  case,  the  mean 
beam  of  reemission  is  at  and  the  particles  are  rccmitted  with  a  Lambertian  distribution.  This  is 
often  characterized  by  saying  that  the  particle  has  lost  all  knowledge  of  its  incoming  direction.  Cook 
[4,13]  fteats  these  two  limiting  cases.  The  actual  function  6X9)  cannot  be  represented  by  such  a  simple 
analytic  function.  The  Hurlbut,  Sherman,  and  NociUa  model  (HSN)  assumes  a  Maxwellian  distribution 
ior  (he  function  'V,  while  leaving  6X9)  unspecified.  Part  of  the  task  of  this  analysis  is  to  determine  this 
function. 
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3.  DRAG  COEFFICIENT’S 


TTie  modem  treatment  of  molecular  scattering  of  gase^  at  a  surface — ^leading  to  definition  of  the 
drag  coefficient  and  lift  coeffideiit  Cp-is  presented  by  Schaaf  lH-  and  [  chaaf  and  Chambre  [15], 
Lacking  detailed  quantum  mechanical  models,  it  is  customary  to  express  the  scattering  in  terms  of  the 
tangential  momentum  transfer  coefficient  c  and  the  ncnnal  momentum  transier  coefficient  a'.  It  begins 
by  assuming  a  gas  in  Maxwellian  equilibrium,  at  temperature  T,  flovring  past  a  surface  element,  with 
velocity  17,  at  an  angle  of  attack  6  (defined  in  Figure  4)  with  the  distribution  function 


(l-u  sine)^*(v*v 


m(2irRr)2 

where  the  gas  constant  R=  k/CN^  M)  =  8.31451  ^xlO’/M  (cm/sec)Vdeg  K,  and  M  is  the  mass  in  atomic 
mass  units  (amu).  Table  1  lists  some  basic  physical  constants  used  in  the  analysis. 


atT414 


Figttne  4.  Coordinate  system  for  surface 
scatterini,  calculations. 
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TABLE  1 

Physical  Constants 


Boltzmann 

Constant 

k  =  1.380658x10'^® 

ergs/K 

Avogadro's 

Number 

N*=  6.0221367x10“ 

mot’ 

Stephan’s 

Constant 

0  =  5.67051x10-® 

W  cm*  K-* 

M(p)  =  1.007276470 

g/mol 

trv  « 1.6726230x10** 

grams/amu 

1.6021773x10’^ 

©rgs/eV 

SchantUwrg  (10, 11]  Introduced  the  Idea  of  cjqjrcssing  the  di-ag  coefficient  in  terms  of  the  thermal 
accommodation  coctfldent  a  Cook  [4],  based  on  Schamberg’s  idea,  provided  the  treatment  most  often 
used  in  aeronomy  studies.  Alfonso  ct  al.  [16]  extended  Schamberg's  ftieory  to  account  for  the  Maxwdlian 
pK^xxties  of  the  gas. 

The  thag  and  lift  ftom  scattering  off  a  flat  surface  element  (Hguie  4)  at  an  incidence  angle  9  can 
be  (^ftainoi  ftom  Sdiaaf  and  Chambre  [IS]  as 


lv/»5  ^^vS 


stn0 


T  S 


-<»)1 


^sinB 

T 
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and 


C,  =  cos^[l+e^(50)|(2-aO(sm^0+^)+acos2e+^^  -^■“ 
+sm0cos0e'(^*)^ 


fl  2V^  a  ^  o' 


\f^S  \firS  2iS^sui0  N  ^ 


wdich  depend  on  o  and  o',  the  speed  ratio  S  =  IJI4' (2RT),  and,  following  Hurlbut  [17],  56  =  5  sin(6). 
Note  ^hat 

^mostpnbaUc^'l’^^  '  ^  ^  pTob^lc  Velocity,  v^=\jW.TIrr  » ^3'  Alfonso  et  al„ 

belcw,  is  the  sveLa%Q  velocity,  and 

is  the  velocity  corresponding  to  the  average  thermal  energy  [18]. 


Scharoberg  [10,11]  ialrodu':d  the  idea  that  Q  could  be  calculated  from  energy  considerations 
aaJ  proposed 

C^=jl+~yr^  sinfljsinO 
with  die  coniesponding  lift  coefficient 

'  3 

Scbamberg's  formulation  assumes  the  condidonr  5»1  and  T„«T.  Schamberg  discusses  a  general 
scattering  pattern  and  adopts  a  Lambertiaa  reemission  model  leading  to  these  relations.  However,  this 
fonnul^on  also  assumes  a  relahon  between  a  and  o  and  o'.  Such  a  relation,  if  indeed  it  exists,  will  be 
more  complex  than  implied  by  these  rslatioDs. 

The  formulas  for  CJ^  for  a  .sphci^  follow,  using 


9 


where  a  is  the  radius  of  the  sphere,  A  =  and  0.  In  this  case  the  integral  in  6  runs  from  -id2 
to  7c/2  to  include  the  interaction  of  thermal  motion  on  die  down  stream  side  of  the  sphere.  The  expression, 
assuming  a  and  o'  are  constant,  is 


(2-o'+o) 


(l+J - L.)erf(Sy(U—) 

4S*  7S^ 


e  ^ 

yfrrS^ 

3SN 

T 

Using  Schamberg's  formulatioo,  we  have 


Here,  the  integral  in  $  runs  from  0  to  Ji/2.  Tlffiy  do  reduce  to  the  same  result  for  the  limiting  case  where 
5»1,  T^«T  and  a  a  o  =  o'  =  1. 

Now  Alfonso  et  al.,  tried  to  include  the  effects  of  the  Maxwell  gas  average  thermal  velocity, 
VfS  / (8RT/X),  in  Sebambestg's  formulation  and  d^ved 


In  the  case  when  a  »  o  =  O  =  1  and  T^«T,  the  Schaaf  and  Chambre  formula  reduces  to  the  slightly 
different  relation 


C**.2*| 


32 


Using  V;,v  for  the  thermal  velodty,  die  relation  becomes 
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,2 


^4 


^Sph_r,4\K 


m>  aurgf 

u 


2j  ^av  energ/ 


Hie  Schaaf  and  Chambre  [IS],  Sdbamberg  [10,11],  and  Alfonso  et  al.  [16]  fonnuladons  cannot  be 
leconciled.  We  adopt  the  Scbaff  and  Chambre  formulation. 


« 


11 


4.  THE  HURLBUT,  SHERMAN,  AND  NOaLLA  THEORY 


In  a  numbo’  of  discussions  of  rare  gas  solid  interactions  for  satellite  regimes  (i.e.,  with  satellite 
velocities  and  tiiamosphere  composition),  Hurlbut  [8,12]  has  summarized  the  existing  data  on  satellite 
accommodation  coefficients.  Nocilla  pointed  out  the  characteristic  of  the  scattered  gas  beii^  a  "drifting 
Maxwellian,"  and  Hurlbut  and  Sherman  [19]  developed  a  formalism  for  malhmnatically  describing  tiiis 
process,  called  the  Hurlbut,  Sherman,  and  Nocilla  (HSN)  tiieory.  In  this  section,  fixe  results  are  given, 
without  s^tproximation,  in  the  form  used  for  the  subsequent  analysis. 


The  HSN  theory  postulates  an  incident  Maxwell  gas  with  velocity  U,  temperature  T,  impinging 
a  surface  at  an  angle  6.  ITie  particles  are  scattered  with  a  velocity  17^  and  ten^rature  at  an  angle 
Given  these  six  variables,  the  molecular  weight,  isenlropic  exponent  y,  and  the  tenqterature  of  the 
scattering  surface  the  HSN  theory  predicts  a,c,  and  o'.  Recall  that  the  adiabatic  exponent  y  =  5.G  for 
point  atoms,  y  =  7/5  for  diatomic  atoms,  and  y  =  4/3  for  polyatomic  atoms  [18].  The  elements  of  tiie 
theory,  as  desoibed  in  Hurlbut  [17],  are  as  follows. 


Following  Scbaaf,  Hurlbut  [17]  finds  the  normal  momentum  flux  p,  of  a  particle  scattered  from 
a  surface 


Pr 


(1-cr  ')|(Sd)e  j(l^Sfl))  j 


'rrT 


2N  T 


This  leads  to  the  expression  for  o'  as 


where 

and 


Now  the  net  tangential  momentum  c  inqxaited  to  scattered  molecules  leads  to 

[T^ 

Tr  «W(®) 
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or 


V  cos(«) 

which  is  somewhat  more  intuitive. 

The  equation  for  a  can  he  obtained  &om  Hurlbut  and  Sherman  [19]  as 


where 

T(S,fl)=(5  ^*^){S9)il^Se)) 

whi(^  for  5»1  reduces  to  \[r(5,6)«>2/tc  ^  (56).  Then  a  is  iq)proximately 


A  useful  result  obtained  when  56»1  and  r^«r  is 


The  resulting  etpuions  provide  a  self-consistent  computation  for  tbs  three  accommodadon 
coefficients  a,  o,  and  o'  in  terms  of  V,  T,  6,  9^  and  y.  These  v^ues  of  the  reflected  parameters 

7%  and  0,  must  be  determined  from  measurements.  Assume  that  there  are  observations  of  O((6}.0(6), 
and  (j'(6)  for  a  given  T,  V,  y,  and  7^  We  can  proceed  as  follows. 

CUmsider  the  three  simultaneous  equations 

[a]f 

[al  aiUJA) 


14 


where  JJ,  T,  6,  and  y  are  known  constants.  Given  the  three  quantities  c,  o',  and  a,  the  values  of 
mi  6, not  independent  but  must  be  determined  as  a  set  A  numerical  inversion  of  these  equations 
can  be  obtained  without  making  the  Si»l  and  5;»1  approximation.  An  inqwrtam  side  result  is  that 
because  and  must  both  be  positive  quantities,  not  aU  combinations  of  o,  o',  and  a  are  possible. 

Start  by  using  the  equation  for  o  to  obtain  £7,  as  a  function  of  6^  and  o,  U  =U(1 

cos  6^ 

Then,  use  this  value  in  the  equation  for  o'  to  obtain  This  is  done  by  iteration.  With  a  little  algebra 
one  can  obtain  the  following  equation,  which  conv^ges  r^idly  by  iteration. 


/ 


JL 


xm-a^se) 


\frr  j~^w 

■  1 


hm 


xm 


U^hid^ 


where 


and 


where  x  contains  the  only  dependence  on  In  effect  we  have  the  function  T,(9,)  withmu  any 
sf^xinutiotL  One  can  now  cono^te  a(£/,(9,),7,(d,),d,},  and  numerically  find  the  value  of  B,  that 
oonesponds  to  the  given  value  of  a  In  doing  this  calculadon  for  the  full  range  of  it  is  immecUately 
ai^ent  that  a  limited  range  of  a  is  obtained  for  a  given  pair  of  or  and  o’. 

As  an  illustration,  we  take  the  example  from  Hurlbut  [17]  for  a  proposed  Magellan  Aerop^  of 
Venus.  Table  2  provides  the  atmosphedc  parameters  for  Venus  at  140-km  altitude. 

Hurlbut  was  illustrating  the  relation  b^een  the  acoommod^on  coefticienis  {a,a,a')  and  the 
scattered  i^rameters  {U^T^BX  His  results  were  based  on  approximations,  where  this  computation  is  done 
withoitt  (qiproximatioa  The  computation  begins  by  dx)0^j3g  the  values  for  ot,  o,  and  d*  given  in  the  first 
three  columns.  The  numerical  inversion  is  performed  as  described  abox-e.  The  first  two  lines  of  the  table 
give  good  inversions.  Added  to  the  table  are  values  of  the  minimum  and  maximum  values  and 
and  the  angle  contending  to  the  latter  6,^).  For  the  other  cases  chasen  by  Hurlbut,  rows  3  through 
8,  the  value  of  a  fa^  outside  the  allowable  range.  Four  additimal  lines  to  die  table  were  added,  choosing 
a  value  of  a  within  the  allowable  range,  and  the  inversion  is  successful. 
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S'  S' 


TABLE  2 
Atmosphoric 
Parameters  for 


This  analy^  leads  to  a  question.  Tbo  range  of  a  could  be  restrided  hiither  by  assuming  that 
T;>T^  Muidvdy.  it  would  seem  plau^ble  that  an  aocomnmdated  molecule  would  always  scatter  with  a 
temperature  ^eater  than  the  scaheiing  surface.  However,  this  seems  to  impose  some  physics  from  what 
fundamentally  a  phenomenological  model.  In  addition  intuitively  oiw  might  believe  that  9;>6.  This 
bourn  out  by  the  data  on  neutral  moleodes  but,  as  will  be  shown  in  Table  3,  not  for  ions. 
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o.es  0.90 


TABLE  3 

HSN  Model  for  Accommodation  Coefficients 


0.90  0.90  0.90 


0.00  0.80 


0.80  0.80  0.90 


0.80  0.80  0.80 


0.90  0.85  0.90 


0.90  0.85  0.95 


0.90  0.85  OJ 


0.90  I  0.90  0.97  S.6 


■ 


S, 

3.185 

0.693 

0.835 

33.7 

1.835 

0.693 

0.835 

33.7 

0.539 

0.968 

0.993 

47.8 

0.381 

0.968 

0.993 

47.8 

1.298 

0.906 

0.959 

39.1 

0.830 

0.906 

0.959 

39.1 

0.630 

0.961 

0.968 

53.7 

0.870 

0.951 

0.988 

63.7 

0.790 

0.959 

0.983 

36.3 

1.118 

0.959 

0.983 

36.3 

0.840 

0.968 

0.993 

47.8 

1.625 

0,906 

0.969 

39.1 

0.825 

0.951 

0.988 

63.7 

1.238 

0.959 

llllllllllll^^ 

36.3 

S.  GOODMAN  AND  WACHMAN  THEORY 


The  Schamberg  theory  for  calculation  of  ballistic  coefficients  uses  the  thermal  accommodation 
coefficient  a  to  diaracterize  the  scattering.  Many  aeronomy  studies  use  the  formulation  by  Cook  [4]  for 
this  calculation.  The  basic  idea  is  as  follows.  From  the  simple  classical  theory,  the  energy  exchange 
between  two  smooth,  hard  spheres  leads  to  the  elementary  resdt 


a- 


4/x 


for  the  mass  ratio  p  =  where  M  is  the  mass  of  the  incident  molecule  and  is  the  mass  of  the 

surface  atom.  Cook  [4]  argues  that  for  most  satellite  surfaces,  tiie  outer  layer  is  primarily  oxygen.  This 
is  either  from  aluminum  oxide  or  silicates  that  are  used  for  solar  cells  or  adsorption  of  atomic  oxygen. 
Therefore,  we  adopt  \6m  the  following  analysis.  In  Baule's  original  analysis  [20],  with  collisions 
between  hard  ^eres,  for  averages  over  all  angles  of  incidence,  the  result  is 


a= 


Goodman  and  Wachman  [21]  derived  a  formula  for  a(7)  based  on  a  lattice  theory  and  experimental  data. 
Though  not  directly  applicable  to  the  satellite  context,  it  gives  a  working  relation  to  be  used  m  the  context 
of  the  Schamberg  formulation.  Goodman  and  Wachman  find  the  relation 


a(oo)® 


for  the  hard  sjdiere  limit  with  p<0.84.  The  final  result  is 


flf(7)=l-«  ^■Hi((oo)tanhl 


a((o)  k 


To 

T 


where  a,  and  X  are  physical  properties  of  the  scattering  surface.  From  Table  II  of  Goodman  and 
Wachman,  we  can  adopt  the  following  values  for  tyi^  satellite  surfaces:  T„=>0,  a  ^  1.4SA‘',  and  X  =: 
670  (g**  deg'Vmolc'^  A).  This  rcthices  to 


a(oo)  k 


For  the  range  of  temperature  considered  here,  we  use  ha  the  following  numerical  examples  for  the 
thermal  accommodation  coefficient 
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6.  EXPERIMENTAL  DATA 


The  Hurlbut,  Sherman,  and  Nodlla  (HSN)  theory  is  based  on  properties  of  the  incident  flux  (C/, 
T,  6)  and  the  reflected  flux  (U,,  T„  6^)  as  well  as  M,  and  y.  and  6,.  would  be  the  underlying 
variables  of  die  scattering  process.  We  would  like  to  have  a  reductionist  model  that  predicts  U,  (U,  T, 
ft  Tw>  Y»  gas  constituents,  scattering  surface  properties),  (U,  T,  9,  y,  gas  constituents,  scattering 
surface  properties),  and  6^(U,  T,  ft  T^,  y,  gas  constituents,  scattering  surface  properties).  Today,  no  such 
model  exists. 

There  are  some  experimental  data  and,  therefore,  we  seek  an  empirical  model  that  is  limited  to 
the  satellite  context  as  defined  above.  The  experimental  data  sets  are  listed  in  Table  4. 


TABLE  4 

Data  Used  in  Model 


Source 

Energy 

(eV) 

Velocfty 

(km/sec) 

o 

c 

Seidel  and  SteinheH 

He 

0.05 

1.77 

-0.001 

0.012 

Musanov  and  Nildforov 

N* 

2.0 

4.2 

0.000 

0.008 

Liu,  Sharma,  Knuth 

He 

1.0 

6.946 

•0.007 

0.046 

Boring  and  Humphris 

N* 

8.0 

7.42 

0.(X)3 

0.021 

Doughty  and 

Schaetzle 

Ar 

25.0 

10.99 

0.001 

0.031 

Boring  and  Humphris 

Nj 

25.0 

13.13 

0.102 

Doughty  and 

Schaetzle 

N* 

25.0 

13.13 

•0.001 

0.026 

Knechtel  and  Pitts 

N,- 

10.0 

8.3 

These  are  all  data  of  (he  momentam  accommodation  coefficients  in  nearly  satellite  conditions.  There  is 
ao  data  of  tbs  thiee  accommodation  coefficsents  together.  We  have  »drxtcd  data  sets  with  velocities 
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corresponding  to  near-earth  orbits.  In  addition,  data  sets  were  limited  to  those  corresponding  to  Cook's 
prescription,  i.e.,  where  a  reasonable  assumption  that  oxygen  was  the  principal  scattering  atom. 

These  data  sets  have  a  number  of  limitations.  First,  there  is  no  data  on  the  thermosphere 
constituents:  N,  O,  H,  and  Oj.  Second,  all  the  results  are  given  in  graphical  form.  Lacking  tabular  data, 
the  numerical  values  of  cr  and  of*  were  necessarily  read  from  the  curves  with  a  consequent  increase  in 
error.  The  numerical  values  are  given  in  Table  11.  The  standard  error  of  fit  of  the  mo^l  to  the  data  is 
0.026,  and  the  mean  and  starulard  error  for  each  data  set  is  given  in  Table  12.  Finally,  the  experimental 
conditions  are  not  thoroughly  known.  For  example,  we  ^  not  know  the  appropriate  values  for  the 
Maxwell  gas  temperature  T  (700  K  was  assumed),  though  we  do  know  the  velocity.  Also,  the  temperature 
of  the  scattering  sample  T*  was  not  given  and  was  assumed  to  be  300  K.  A  brief  description  of  each  data 
set  follows. 

6a  SEIDEL  AND  STEINHHL  (1974) 

The  description  of  the  experimental  set  up  stated  that  the  He  velocity  was  1770  m/sec.  However, 
it  was  also  stated  that  tWs  corresponds  to  an  energy  of  0.05  eV.  He  at  1 .770  km/sec  has  a  energy  of  0.065 
eV,  and  He  with  an  energy  of  0.05  eV  has  a  velocity  of  1^53  km/sec.  We  have  chosen  to  characterize 
the  data  with  a  velocity  of  1.770  km/sec. 

Seidel  and  Steinheil  chose  to  give  their  results  in  terms  of  a  and  the  ‘relative  nonnai  momentum 
transfer  RNT“  ct,'.  For  use  here,  it  was  converted  to  o"  =  2-€t,'‘. 

Seidel  and  Steinbdl  piotted  results  fw  He  scattering  off  copper,  shellac,  tong^n.  gold,  gto.  and 
sapf^re.  FQQowing  Cook's  reasoning,  we  sel»;ted  the  saj^ihire  data  as  a  sapphke.  AljO),  surface  is  mo^ 
Ul^y  to  present  oxygen  as  Us  ^Meting  object 

These  data  are  used  to  obtain  an  initial  evaluation  of  tl»  reflected  scatting  angle  Using  tire 
approximate  formula  for  0„  we  can  calculate  6,  from  o.cf.  and  9.  As  seen  in  Figure  5.  for  He  at  1.77 
km/$ec.  6,  is  a  monotooically  increasing  friiK:tion  of  9.  At  incident  angles  ;^oacMi^  90  deg.  it  is  nearly 
equal  to  9.  This  is  the  condition  tor  specular  reflection. 

6,2  MUSANOV,  NIKIFOROV,  OMEUK,  AND  FREEDLENDER  (19$5) 

These  authors  (22)  provide  data  for  Nj  scattering  off  of  aluminum,  duralium,  steel,  and  a  variety 
of  tin  surfaces.  The  Nj  velocity  is  given  as  2  cV  or  4.7  km/scc.  Tivey  claim  a  vcloctiy  ratio  of  S  = 
5-5iD.7,  whidi  would  suggest  a  Maxwell  tcmpcratuire  of  7  =  982  K.  Because  of  the  uncertainty  of  S,  we 
have  adopted  T  ^  7()0  K  consistent  with  the  other  data  sets.  The  autbcNrs  also  state  that  T,  3IQ±20  K. 
Because  of  the  uncertainty,  we  have  adc^tted  T.  s  3(X)  K  consistent  with  the  otha  data  sets. 


J31741-5 
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Figure  5.  Scattering  angle:  helium  on  sapphire. 


Musanov  et  al.  have  reduced  their  data  and  provide  an  analytical  expression  that  can  be  used 
to  obtain  c  and  a\  We  have  chosen  the  data  for  aluminum,  assuming  that  it  is  in  fact  AlOj.  The 
eiqiressions  are 

<j  =  1.01 


and 


o'  =  0.5906-0.1044  cos(2e)  . 

'These  expressions  were  evaluated  between  0  aiuJ  90  deg  at  10-deg  intervals.  These  become  the  tabular 
dat'^  used  in  the  least  squares  determination  of  model  parameters. 


63  LIU,  SHARMA,  .\ND  KNUTH  (1979) 

These,  exjieriments  used  1  eV  (6,946  kro/sec)  He  particles  imp'uging  on  an  aluminum  plate  and 
an  anodized  elummum  plate.  They  stuthed  scattering  as  a  t>mction  of  incid^t  angle  in  terms  of  the 
dis.tibution  of  scattered  particles  in  the  plane  containing  the  incident  beam  and  the  surface  normal  as  well 
as  the  distribution  out  of  this  plane.  They  adopted  the  representation  of  etjv  =  to  represent  the 

data.  The  latter  were  converted  to  &  using  “f,,,  =  300  K. 

6A  KNUTH  (19(10) 

Kauth  did  not  publish  any  new  raeasureuwnts.  However,  he  converted  a  number  of  existing 
normal  momentum  accocimodaho''-  coefficient’s  measuremoa,;  into  the  variable  a"  =  and  published 
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th  .!  numerical  values.  The  purpose  was  to  demonstrate  that  a  single  model  could  be  used  for  this  variable, 
i,e.,  a"  =  1  +  0/90, 6  expressed  in  degrees.  He  did  provide  a  table  of  values  talren  from  the  same  authors 
desaibed  here.  With  the  view  of  examining  this  simple  relation,  this  data  were  also  included  in  the 
determination  of  the  HSN  model  parameters. 

6S  BORING  AND  HUMPHRIS  (19^.^0) 

This  experimental  data  [23]  wns  of  in  the  energy  range  8-200  eV,  7-37  km/sec.  The 
aq)eiime.nts  were  done  on  samples  of  the  material  used  for  the  ^ho  I  and  Echo  n  satellites.  Echo  I  was 
alumimun  evaporated  on  Myiar,  and  Echo  n  was  aluminum  with  a  coating  of  Alodine,  an  amorphous 
phosphate  The  Echo  I  data  was  choser?  ibr  this  aitalysis.  Data  was  read  from  the  curves  at  8  eV,  7.4 
km/sec  and  25  eV,  13.13  km/sec,  die  higher  velocities  being  much  greater  than  an  expected  satellite 
velocity. 

Boring  and  Humphris  provide  the  ohserv-able 


p 

-/^=(2-a-o^cos^6-(l-c)+0.055  cos0 


This  is  used  directly  in  the  confutation  of  the  HSN  model  parameters. 


6,6  DOUGHTY  AND  SCHAETZLE  (1969) 

Doughty  and  Schactzle  [24]  made  measurements  of  Nj,  air,  and  At  covering  the  energy  range 
4-200  eV.  Sdiactide  [25]  describes  the  cxperimcmal  setup.  Their  scattering  materials  included  aluminum 
and  fresh  varnish.  Only  data  on  normal  scattering,  0  =  90  deg,  was  given  covering  die  full  energy  range. 
The  data  for  different  incidence  angles  were  given  for  25, 75,  and  150  eV.  Only  the  25  cV  data  (10.99 
km/sec  for  Ar  and  13.13  km/sec  for  Nj)  are  used  in  this  analysis.  All  of  this  data  exhibits  the 
characteristic  o  1.0  for  some  range  of  9,  This  leads  to  "frirward  scattering,"  l.e.,  0,>9O  deg  where  die 
sca^ered  molecule  reverses  direction.  This  is  shown  in  Eigure  6.  Note  the  range  fw  0,  is  0  to  180  deg 
and  that  0,  =  112.6  deg  0  =  60  deg  for  both  Nj  and  Ar,  In  other  words,  the  lobe  of  the  reflexion 
pauem  is  bent  rearward  rather  than  in  the  forward  dire<^on. 
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Figure  6.  Scattering  angle:  aluminum. 


6.7  KNECHTEL  AND  PITTS  (1969,  1973) 

The  data  [26,27]  are  taken  with  Ni*  at  5, 10, 15,  and  20  eV  (5.87, 8.30, 10.17,  and  11.74  km/sec). 
The  scattering  surface  is  aiuminmn.  Figure  7  plots  the  reflected  scattering  angle  6^  The  behavior  of  6, 
is  quite  different  than  the  two  previous  examples.  We  found  that  the  ion  data  could  not  be  made 
consistent  with  the  other  data  used  in  this  analysis,  and  the  data  were  not  included  in  the  determination 
of  the  HSN  model  parameters. 


figure  7.  Scattering  angle:  Nj*  on  aluminum. 
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7.  DATA  ANALYSIS 


The  data  analysis  began  by  computing,  where  possible,  6^  from  the  approximate  formula.  The 
evidenced  backscattering  mentioned  by  Hurlbut  is  seen.  This  is  particularly  evident  in  the  high-velocity 
measurements  of  Doughty  and  Schaetde  where  tiie  scattering  an^e  exceeds  90  deg.  In  general,  for  each 
data  set,  6,  inaeases  monotonically  with  0  and  0p»9.  However,  the  data  for  N/  is  anomalous  in  diat  6,<Q 
and  monotonically  deaeases.  The  data  on  N2+  were  not  used  in  the  analysis.  We  assume  that  this  model 
is  inappropriate  for  ions. 

Motivated  by  the  Baule  model  for  ot,  an  attempt  was  made  to  find  a  relation  in  terms  of  the  mass 
ratio.  In  this  data,  the  mass  ranges  from  4  to  40  amu.  However,  the  large  difference  in  mass  for  Nj  and 
Ar,  and  the  small  diffidence  in  o  and  o'  for  the  two  constituents,  frustrated  this  model  dependence. 

RnaUy,  a  model  that  depended  only  U,  0,  and  the  mass  ratio  p  was  found  to  be  acceptable.  The 
model  dependence  adopted  is  then  -  U/U,6,n),  T,  =  T^U,6),  and  =  6/U,6).  This  was  done  by 
referring  all  the  constituents  to  He.  The  He  model  was  in^jlemented  as  follows.  A  grid  of  points  (UpBj) 
is  defined  separately  for  U„  T^,  and  B^  For  a  given  f/,0  combination,  linear  interpolation  is  used  within 
this  grid,  i.e.,  a  ruled  surface.  The  values  of  the  functions  ({7„  and  B^  at  the  grid  points  are  determined 
by  a  least  squares  conq)utation,  as  described  below.  Then,  the  values  of  17,  for  die  constituent  x  is 
computed  from 


(  ir*  \- 

JZ-  2_1 

..Ht 
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where  p  is  conqiuted  assuming  oxygen  is  the  scattering  atom.  The  values  for  T,  and  B,  are  used  without 
correctioiL 

The  grid  point  values  were  determined  widi  an  iterative  least  squares  program  that  accepts  as  input 
the  grid  point  model  and  observatioits  of  ct,  o,  o', a",  and  the  linear  combination  observed  by  Boring  and 
Humphris  (2-o-o')cos^0-(l-o)+O.O55cos0.  This  nonlinear  least  squares  conqiutation  was  iterated  until 
convergence.  It  is  found  that,  because  of  the  large  value  of  the  speed  ratio  S  =  UI-fi2RT),  the  coefficients 
for  1/  and  r  are  almost  conqiletely  correlated  numerically.  This  was  overcome  by  tqiplying  weak  prior 
knowledge  [28],  using  an  initial  variance  of  the  temperature  variable  of  ±1000  K.  This  required  choosing 
the  initial  values  of  T,  which  was  done  by  experiment.  In  addition,  from  symmetry  considerations,  the 
values  of  BX^  =  90  deg)  were  constrained  to  be  90.0  deg.  The  final  least  squares  solution  used  105 
observations,  had  a  mean  error  of  -0.003,  and  a  standard  error  of  0.02S.  The  total  data  set  fit  this  model 
with  an  accuracy  of  2.5%.  Table  5  lists  the  grid  values  for  U„  Table  6  lists  the  grid  values  for  Tp  and 
Table  7  lists  the  grid  values  for  6^  Tables  8,  9,  and  10  give  the  formal  standard  error  a  for  each  grid 
point 
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TABLE  5 

Grid  Values  for  (km/sec) 


U(km/sec)\e 

0.0 

45.0 

90.0 

M7 

2.146 

0.840 

0.225 

4.20 

2.806 

4.041 

3.989 

6.946 

4.020 

3.143 

3.174 

10.990 

7.367 

7.557 

5.649 

13.130 

7.744 

4.944 

4.474 

TABLE  8 

Grid  Values  for  7,  (K) 


U(km/sec}\9 

0.0 

45.0 

90.0 

1.77 

1 

35 

75 

13.130 

992 

1003 

997 
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TABLE  7 

Grid  Values  for  Sf  (deg) 


U(km/sec)\0 

0.0 

30.0 

45.0 

60.0 

75.0 

90.0 

t.770 

74.3 

70.1 

65.8 

67.1 

71.4 

90.0 

4.200 

96.5 

91.2 

91.3 

90.8 

90.4 

90.0 

6.946 

89.4 

91.6 

90.8 

90.6 

89.7 

90.0 

7.400 

53.8 

61.1 

62.4 

54.3 

62.3 

90.0 

10.990 

■ 

65.0 

104.1 

111.3 

128.4 

90.0 

13.130 

0.0 

67.3 

81.6 

125.8 

127.5 

90.0 

TABLES 


Grid  Value  a  for  (km/sec) 


U(kiTVsec)\6 

0.0 

45.0 

90.0 

1.77 

0.09 

0.08 

0.13 

4.20 

0.02 

0.04 

0.07 

6.946 

0.11 

0.08 

0.08 

10.990 

0.49 

0.10 

0.16 

13.130 

llllllllQl 

0.12 

0.18 

TABLE  9 

Grid  Valutt  o  for  T,  (K) 


U(km/sec)\e 

IH 

45.0 

90.0 

1.77 

24 

■ 

22 

13.130 

25 

25 

25 

TABLE  10 

Grid  Vatuo  a  for  0^  (dag) 


U(krT^sec)\6 

■ 

30.0 

45.0 

60.0 

75.0 

90.0 

1.770 

5.6 

1.8 

m 

3.1 

6.2 

0.0 

4.200 

10.6 

4.1 

3.1 

1.8 

■ 

lllllllllllll 

6.946 

m 

2.5 

2.2 

1.6 

0.8 

m 

7.400 

10.8 

6.3 

9.8 

8.1 

12.2 

Ifl 

10.990 

27.6 

3.8 

3.8 

■ 

9.5 

0.0 

13.130 

24.3 

4.3 

5.5 

4.1 

10.6 

o 

b 

The  M  data  set  and  residuals  ate  given  in  Table  11.  The  re^duals  for  the  Knechtel  and  Pitts 
data  are  also  Induded,  though  they  were  not  used  in  the  least  squares  adjustment  It  does  not  fit  with  this 
model.  Table  12  provides  the  statistics  for  each  dsda  set  sepaikely. 
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Sp  U  T 

km/s  K 


TABLE  11 

Residuals  for  Least  Squares  Fit 


Tw  c  th  obs 
K 


N2 

4.200 

700 

300 

n 

90 

0.695 

N2 

4.200 

700 

300 

t 

90 

1.010 

N2 

4.200 

700 

300 

n 

80 

0.689 

N2 

4.200 

700 

300 

t 

80 

1.010 

N2 

4.200 

700 

300 

n 

70 

0.671 

N2 

4.200 

700 

300 

t 

70 

1.010 

N2 

4.200 

700 

300 

n 

60 

0.643 

N2 

4.200 

700 

300 

t 

60 

1.010 

N2 

4.200 

700 

300 

n 

50 

0.609 

N2 

4.200 

700 

300 

t 

50 

1.010 

N2 

4.200 

700 

300 

n 

40 

0.572 

N2 

4.200 

700 

300 

t 

40 

1.010 

N2 

4.200 

700 

300 

n 

30 

0.538 

N2 

4.200 

700 

300 

t 

30 

1.010 

N2 

4.200 

700 

300 

n 

20 

0.511 

N2 

4.200 

700 

300 

t 

20 

1.010 

N2 

4.200 

700 

300 

n 

10 

0.492 

N2 

4.200 

700 

300 

t 

10 

1.010 

N2 

4.200 

700 

300 

n 

0 

0.486 

N2 

4.200 

700 

300 

t 

0 

1.010 

N2 

7.420 

700 

300 

b 

90 

0.160 

N2 

7.420 

700 

300 

b 

75 

C.130 

N2 

7.420 

700 

300 

b 

60 

0.100 

N2 

7.420 

700 

300 

b 

30 

0.040 

N2 

7.420 

700 

300 

b 

15 

-0.020 

N2 

7.420 

700 

300 

b 

90 

0.195 

N2 

7.420 

700 

300 

b 

75 

0.137 

N2 

7.420 

700 

300 

b 

60 

0.080 

N2 

7.420 

700 

300 

b 

45 

0.070 

N2 

7.420 

700 

300 

b 

30 

0.001 

N2 

13.130 

700 

300 

b 

60 

0.070 

N2 

13.130 

700 

300 

b 

30 

0.030 

N2 

13.130 

700 

300 

b 

15 

-0.040 

N2 

13.130 

700 

300 

b 

90 

0.100 

N2 

13,130 

700 

300 

b 

50 

0.050 

N2 

13.130 

700 

300 

b 

45 

0.020 

N2 

13.130 

700 

300 

b 

30 

-0.010 

N2 

13.130 

700 

300 

b 

15 

-0.130 

N2 

13.130 

700 

300 

n 

90 

0,857 

N2 

13.130 

700 

300 

n 

90 

0.886 

N2 

13.130 

700 

300 

m 

90 

1.140 

res 

o 

a’ 

a" 

a 

-0.003 

0.646 

0.698 

1.363 

0.879 

0.001 

1.009 

0.698 

1.363 

0.879 

-0.002 

1.010 

0.691 

1.371 

0.878 

0.000 

1.010 

0.691 

1.371 

0.878 

-0.001 

1.010 

0.672 

1.391 

0.877 

0.000 

1.010 

0.672 

1.391 

0.877 

0.004 

1.010 

0.639 

1.426 

0.876 

-0.000 

1.010 

0.639 

1.426 

0.876 

0.024 

1.01  i 

0.585 

1.482 

0.874 

-0.001 

1.011 

0.585 

1.482 

0.874 

0.021 

1.010 

0.551 

1.524 

0.896 

0.000 

1.010 

0.551 

1.524 

0.896 

-0.015 

1.007 

0.553 

1.541 

0.934 

0.003 

1.007 

0.553 

1.541 

0.934 

-0.026 

1.011 

0.537 

1.590 

0.964 

-0.001 

1.011 

0.537 

1.590 

0.964 

0.001 

1.012 

0.491 

1.694 

0.986 

4).002 

1.012 

0.491 

1.694 

0.986 

0.030 

1.008 

0.456 

1.843 

1.000 

0.002 

1.008 

0.456 

1.843 

1.000 

-0.001 

0.763 

0.894 

1.151 

0.973 

-0.004 

0.754 

0.896 

1.151 

0.971 

0.010 

0.829 

0.887 

1.165 

0.969 

0.020 

0.894 

0.710 

1.361 

0.955 

0.037 

0.875 

0.316 

1.742 

0.942 

0.034 

0.763 

0.894 

1.151 

0.973 

0.004 

0.754 

0.896 

1.151 

0.971 

-0.010 

0.829 

0.887 

1.165 

0.969 

0.000 

0.898 

0.836 

1.223 

0.966 

-0.019 

0.894 

0.710 

1.361 

0.955 

-0.140 

1.204 

0.851 

1.177 

0.961X 

-0.026 

0.885 

0.544 

1.487 

0.925X 

0.163 

0.717 

0.307 

1.726 

0.884X 

-0.092 

1.376 

0.863 

1.162 

0.969X 

-0.160 

1.204 

0.851 

1.177 

0.961X 

-0.121 

0.962 

0.757 

1.273 

0.957X 

-0,066 

0.885 

0.544 

1.487 

0.925X 

0.073 

0.717 

0.307 

1.726 

0.884X 

-0.006 

0.850 

0.863 

1.162 

0.969 

0.023 

0.850 

0.863 

1.162 

0.969 

-0.022 

0.850 

0.863 

1.162 

0.969 

31 


TABLE  11  (Continued) 
Residuals  for  Least  Squares  Fit 


Sp 

U 

T 

Tw 

c 

di 

obs 

res 

a 

& 

a" 

a 

km/s 

K 

K 

N2 

13.130 

700 

300 

n 

80 

0.857 

-0.010 

1.385 

0.867 

1.158 

0.966 

N2 

13.130 

700 

300 

n 

75 

0.887 

0.008 

1.382 

0.879 

1.147 

0.965 

N2 

13.130 

700 

300 

m 

75 

1.140 

-0.007 

1.382 

0.879 

1.147 

0.965 

N2 

13.130 

700 

300 

n 

70 

0.857 

-0.014 

1.292 

0.871 

1.155 

0.964 

N2 

13.130 

700 

300 

n 

60 

0.824 

-0.027 

1.204 

0.851 

1.177 

0.961 

N2 

13.130 

700 

300 

n 

60 

0.838 

-0.013 

1.204 

0.851 

1.177 

0.961 

N2 

13.130 

700 

300 

m 

60 

1.190 

0.013 

1.204 

0.851 

1.177 

0.961 

N2 

13.130 

700 

300 

t 

60 

1.237 

0.033 

1.204 

0.851 

1.177 

0.961 

N2 

13.130 

700 

300 

n 

50 

0.764 

-0.017 

1.031 

0.781 

1.248 

0.958 

N2 

13.130 

700 

300 

t 

50 

1.000 

-0.031 

1.031 

0.781 

1.248 

0.958 

N2 

13.130 

700 

300 

n 

45 

0.792 

0.035 

0.962 

0.757 

1.273 

0.957 

N2 

13.130 

700 

300 

m 

45 

1.240 

-0.033 

0.962 

0.757 

1.273 

0.957 

N2 

13.130 

700 

300 

n 

40 

0.676 

-0.025 

0.938 

0.701 

1.330 

0.947 

N2 

13.130 

700 

300 

t 

40 

0.979 

0.041 

0.938 

0.701 

1.330 

0.947 

N2 

13.130 

700 

300 

n 

30 

0.538 

-0.006 

0.885 

0.544 

1.487 

0.925 

N2 

13.130 

700 

300 

n 

30 

0.551 

0.007 

0.885 

0.544 

1.487 

0.925 

N2 

13.130 

700 

300 

m 

30 

1.480 

-0.007 

0.885 

0.544 

1.487 

0.925 

N2 

13.130 

700 

300 

t 

30 

0.875 

-0.010 

0.885 

0.544 

1.487 

0.925 

N2 

13.130 

700 

300 

t 

20 

0.781 

0.011 

0.770 

0.387 

1.645 

0.899 

Hs 

1.770 

700 

300 

m 

90 

1.260 

0.020 

0.873 

1.253 

1.240 

1.223 

He 

1.770 

700 

300 

m 

80 

1.250 

•0.009 

0.747 

1.226 

1.259 

1.220 

He 

1.770 

700 

300 

t 

80 

0.740 

-0.007 

0.747 

1.226 

1.259 

1.220 

He 

1.770 

700 

300 

m 

70 

1.270 

-0.014 

0.718 

1.204 

1.284 

1.215 

He 

1.770 

700 

300 

t 

70 

0.730 

0.012 

0.718 

1.204 

1.284 

1.215 

He 

1.770 

700 

300 

m 

60 

1.300 

-0.019 

0.721 

1.176 

1.319 

1.206 

He 

1.770 

700 

300 

t 

60 

0.720 

-0.001 

0.721 

1.176 

1.319 

1.206 

He 

1.770 

700 

300 

m 

50 

1.380 

0.013 

0.727 

1.135 

1.367 

1.194 

He 

1.770 

700 

300 

t 

50 

0.720 

-0.007 

0.727 

1.135 

1.367 

1.194 

He 

1.770 

700 

300 

m 

40 

1.480 

0.014 

0.719 

1.009 

1.466 

1.159 

He 

1.770 

700 

300 

t 

40 

0,720 

0.001 

0.719 

1.009 

1.466 

1.159 

He 

1.770 

700 

300 

m 

30 

1.620 

-0.015 

0.717 

0.739 

1.635 

1.084 

He 

1.770 

700 

300 

t 

30 

0.720 

0.003 

0.717 

0.739 

1.635 

1.084 

He 

1.770 

700 

300 

m 

20 

1.840 

0.002 

0.701 

0.361 

1.838 

0.982 

He 

1.770 

700 

300 

t 

20 

0,700 

•0.001 

0.701 

0.361 

1.838 

0.982 

He 

6.946 

700 

296 

a 

90 

0.618 

0.026 

0,543 

0.592 

1.489 

0.759 

He 

6.946 

700 

296 

m 

90 

1.470 

-0.019 

0.543 

0.592 

1.489 

0.759 

He 

6.946 

700 

296 

n 

75 

0.574 

-0.005 

0.991 

0.579 

1.503 

0.762 

He 

6.946 

700 

296 

m 

75 

1.510 

0.007 

0.991 

0.579 

1.503 

0.762 

He 

6.946 

700 

296 

t 

75 

0.990 

-0.001 

0.991 

0.579 

1.503 

0.762 

He 

6.946 

700 

296 

n 

60 

0.514 

•0.014 

1.010 

0.528 

'  1.555 

0.765 

He 

6.946 

700 

296 

m 

60 

1.570 

0.015 

1.010 

0.528 

1.555 

0.765 

He 

6.946 

700 

296 

t 

60 

1.010 

0.000 

1.010 

0.528 

1.555 

0.765 
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TABLE  11  (Continued) 
Residuals  for  Least  Squares  Fit 


Sp  u 

T 

Tw 

c 

th 

obs 

res 

a 

a' 

a" 

a 

km/s 

K 

K 

He  6.946 

700 

296 

n 

45 

0.375 

-0.040 

1.009 

0.415 

1.664 

0.767 

He  6.946 

700 

296 

m 

45 

1.700 

0.036 

1.009 

0.415 

1.664 

0.767 

He  6.946 

700 

296 

t 

45 

1.010 

0.001 

1.009 

0.415 

1.664 

0.767 

He  6.946 

700 

296 

m 

30 

1.800 

-0.155 

1.016 

0.060 

1.955 

0.732 

He  6.946 

700 

296 

t 

30 

1.010 

-0.006 

1.016 

0.060 

1.955 

0.732 

He  6.946 

700 

296 

m 

15 

2.620 

0.049 

1.00 

-0.932 

2.571 

0.693 

He  6.946 

700 

296 

t 

15 

1.000 

0.000 

1.000 

-0.932 

2.571 

0.693 

At  10.990 

700 

300 

n 

90 

0.896 

-0.001 

0.884 

0.897 

1.128 

0.980 

At  10.990 

700 

300 

m 

90 

1.100 

-0.028 

0.884 

0.897 

1.128 

0.980 

At  10.990 

700 

300 

n 

90 

0.926 

0.029 

0.884 

0.897 

1.128 

0.980 

Ar  10.990 

700 

300 

n 

80 

0.863 

-0.010 

1.386 

0.873 

1.152 

0.969 

At  10.990 

700 

300 

n 

75 

0.887 

0.013 

1.418 

0.874 

1.152 

0.963 

At  10.990 

700 

300 

m 

75 

1.140 

-0.012 

1.418 

0.874 

1.152 

0.963 

Ar  10.990 

700 

300 

n 

70 

0.813 

-0.030 

1.306 

0.843 

1.182 

0.956 

At  10.990 

700 

300 

n 

60 

0.736 

-0.032 

1.169 

0.768 

1.257 

0.940 

At  10.990 

700 

300 

n 

60 

0.745 

-0.023 

1.169 

0.768 

1.257 

0.940 

At  10.990 

700 

300 

m 

60 

1.280 

0.023 

1.169 

0.768 

1.257 

0.940 

At  10.990 

700 

300 

t 

60 

1.193 

0.024 

1.169 

0.768 

1.257 

0.940 

At  10.990 

700 

300 

n 

50 

0.654 

-0.026 

1.120 

0.680 

1.345 

0.920 

Ar  10.990 

700 

300 

t 

50 

1.111 

-0.009 

1.120 

0.680 

1.345 

0.920 

Ar  10.990 

700 

300 

n 

45 

0.656 

0.035 

1.100 

0.631 

1.404 

0.910 

Ar  10.990 

700 

300 

m 

45 

1.370 

-0.034 

1.100 

0.621 

1.404 

0.910 

Ar  10.990 

700 

300 

a 

40 

0.577 

0.006 

1.007 

0.571 

1.454 

0.911 

Ar  10.990 

700 

300 

t 

40 

0.996 

-0.011 

1.007 

0.571 

1.454 

0.911 

Ar  10.990 

700 

300 

a 

30 

0.495 

-0.010 

0.862 

0.505 

1.524 

0.913 

Ar  10.990 

700 

300 

a 

30 

0.509 

0.004 

0.862 

0.505 

1.524 

0.913 

Ar  10.990 

700 

300 

ai 

30 

1.520 

-0.004 

0.862 

0.505 

1.524 

0.913 

Ar  10.990 

700 

300 

t 

30 

0.864 

0.002 

0.862 

0.505 

1.524 

0.913 

Ar  10.990 

700 

300 

t 

20 

0.710 

-0.073 

0.783 

0.459 

1.579 

0.915 

N2+  8.300 

700 

300 

a 

90 

0.902 

0.059 

0.819 

0.843 

1.195 

0.957X 

N2+  8.300 

700 

300 

a 

80 

0.902 

0.076 

0.853 

0.826 

1.212 

0.95  IX 

N2+  8.300 

700 

300 

a 

75 

0.896 

0.080 

0.846 

0.816 

1.222 

0.948X 

N2+  8.300 

700 

300 

a 

70 

0.884 

0.079 

0.843 

0.805 

1.233 

0.945X 

N2+  8.300 

700 

300 

a 

60 

0.854 

0.074 

0.832 

0.780 

1.261 

0.937X 

N2+  8.300 

700 

300 

t 

60 

0.427 

-0.405 

0.832 

0.780 

1.261 

0.937X 

N2+  8.300 

700 

300 

a 

50 

0.787 

0.063 

0.878 

0.724 

1.319 

0.930X 

N2+  8.300 

700 

300 

t 

50 

0.506 

-0.372 

0.878 

0.724 

1.319 

0.930X 

N2+  8.300 

700 

300 

a 

40 

0.701 

0.051 

0.879 

0.650 

1.395 

0.922X 

N2+  8.300 

700 

300 

t 

40 

0.530 

•0.349 

0.879 

0.650 

1.395 

0.922X 

N2+  8.300 

700 

300 

a 

30 

0.585 

0.040 

0.851 

0.545 

1.503 

0.915X 

N2+  8.300 

700 

300 

t 

30 

0.494 

-0.357 

0.851 

0.345 

1.503 

0.915X 

N2+10.165 

700 

300 

a 

90 

0.970 

0.210 

0.740 

0.760 

1.268 

0.923X 
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TABLE  11  (Continued) 
Residuals  for  Least  Squares  Fit 


Sp  u 

T 

Tw 

c 

th 

obs 

res 

a 

o' 

o" 

a 

km/s 

K 

K 

N2+10.165 

700 

300 

n 

SO 

0.970 

0.237 

1.450 

Q.m 

1.294 

0.905X 

N2+10.165 

700 

300 

n 

75 

0.963 

0.237 

1.472 

0.726 

1.302 

0.895X 

N2+10.165 

700 

300 

t 

75 

0.232 

-1.240 

1.472 

0.726 

1.302 

0.895X 

N2+10.165 

700 

300 

n 

70 

0.963 

0.272 

1.298 

0.691 

1.336 

0.885X 

N2+10.165 

700 

300 

t 

70 

0.378 

-0.920 

1.298 

0.691 

1.336 

0.885X 

N2+10.165 

700 

300 

n 

60 

0.933 

0.324 

1.103 

0.609 

1.417 

0.862X 

N2+10.165 

700 

300 

t 

60 

0.585 

-0.518 

1.103 

0.609 

1.417 

0.862X 

N2+10.165 

700 

300 

n 

50 

0.896 

0.387 

1.061 

0.509 

1.515 

0.838X 

N2+10.165 

700 

300 

t 

50 

0.677 

-0.384 

1.061 

0.509 

1.515 

0.838X 

N2+10.165 

700 

300 

n 

40 

0.841 

0.458 

0.948 

0.383 

1.639 

0.824X 

N2+10.165 

700 

300 

t 

40 

0.689 

-0.259 

0.948 

0.383 

1.639 

0.824X 

N2+i0.165 

700 

300 

n 

30 

0.738 

0.461 

0.794 

0.277 

1.743 

0.824X 

N2+10.165 

700 

300 

t 

30 

0.683 

-0.111 

0.794 

0.277 

1.743 

0.824X 
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In  Table  11  fte  column  labeled  c  indicates  the  observation  type.  An  "n"  denotes  an  observation  of  o',  a 
"t"  denotes  an  observation  of  o,  an  "m"  denotes  an  observation  of  o",  and  a  "b"  indicates  an  observation 
of  the  linear  combination  of  o  and  o'  given  by  Boring  and  Humphris.  The  calculated  values  of  o,  o',  o", 
and  a  are  also  given.  An  "x"  at  the  end  of  the  line  indicates  that  this  observation  was  not  used  in  the 
determination  of  the  model  parameters. 


TABLE  12 

Individual  Statistics  for  Solution 


Specie  Velocity  <mean> 


-0.001 


-0.051 


-0.003 


Author 


0.012  Seidel  and  Steinheil 


0.008  Musanov  and 
Nikiforov 


0.046  Liu,  Sharma,  and 
Knuth 


0.021  Boring  and  Humphris 


0.031  Doughty  and 
Schaetde 


0.026  Doughty  and 
Scha^e 


0.102  Boring  and  Humphris 


0.025  Ail 


As  an  aid  to  viewing  the  scattering  results.  Figures  8,  9,  and  10  show  the  two-dimensional 
dqtendence  of  and  0,.  on  C/  and  0.  The  scattering  angle  surface  reveals  three  regions.  For 
inbsimediate  6>U>^  km/sec,  0,  is  between  60  and  90  deg,  approaching  a  classical  diffuse  scattering  regime. 
For  lower  velocities,  it  tends  toward  specular  scattolng.  For  higher  velocities  the  scattering  exhibits 
forward  scattering  where  0;>9O  deg.  Also  for  higher  velocity,  and  6<4S  deg,  foe  scattering  angle  seetits 
to  follow  0,e20.  The  scattering  velocity  exhibits  a  different  taxonomy.  The  scattering  velocity,  at  grazing 
incidence  0=0,  is  well  correlated  but  less  than  the  incident  velocity. 


At  6=90  deg,  the  scattered  velocity  C/,  has  a  maximum  at  about  17=4  km/sec,  a  minimum  at  U=8 
km/sec,  and  a  maTimiim  again  at  17=11  km/sec.  At  U=8  km/sec  the  velocity  is  relatively  constant, 
which  is  consistent  widi  diffuse  scattering. 


231741<6 


Figurg  <5.  Scaaemg  angle:  helium  on  oxygen. 
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SCATTER.NG  TEMPERATURE  (K) 


23174M0 


Figure  10.  Scattering  temperature:  helium  on  oxygen. 
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8.  CALCULATION  OF  Q 


8.1  COMPASISON  OF  GOODMAN  AND  WACHMAN  THEORY  AND  HSN  THEORY 

The  two  formalisms  for  computing  a  can  hardly  be  erqrected  to  give  the  same  result  as  they  are 
based  on  q>iite  cfifferent  variables.  The  thermal  accommodation  coefficient  of  Goodman  and  Wactoan 
(G&W)  depends  on  the  mass  ratio  in  the  ten^rature  regime  considered  here.  The  HSN  thermal 

accommodation  coefficient  depends  primarily  on  the  incident  velocity  U  and  secondarily  on  p,  T,  and 
y.  The  HSN  theory  also  includes  the  effects  of  a  molecule's  internal  energy,  which  is  assumed  to  be 
undmnged  in  the  accommodating  process.  In  addition,  G&W  give  the  range  of  validity  for  the  theory  as 
p<0.83.  In  Figure  11  we  plot  a  derived  from  the  two  formulations.  The  HSN  values  are  obtained  with 
0=9C  deg.  The  figure  uses  the  mass  ratio  p  as  the  independent  variable  and  plots  the  values  of  a.  The 
plotted  points  are  in  order  of  increasing  p;  H,  He,  N,  O,  N^,  O2,  and  Ar.  The  G&W  values  are 
consistently  smaller  than  the  HSN  values.  There  is  an  observation  of  Oq  by  Krech  et  al.  [29].  They 
measured  a  by  scattering  atomic  oxygen  off  of  three  substances:  RCG,  nickel,  and  gold.  It  is  argued  that 
the  atomic  oxygen  beam  rapidly  deposits  an  adsorbed  layer  of  0  atoms  on  the  surface,  and  that  the 
scattermg  surface  is  then  effectively  O;  i.e.,  p=l.  Krech  et  al.  give  a  value  of  (»=0.64±0.1,  which  is  in 
close  agreement  with  the  value  firom  the  G&W  theory. 


Figure  11.  Comparison  of  G&W  with  HSN. 
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CALCULATION  OF  C4  FOR  FLAT  PLATES 


Uie  HSN  model  obtained  in  Section  7  can  now  be  used  for  calculation  of  Q.  First,  exan:q)les  are 
given  for  flat  plates.  In  Hgure  12  data  are  plotted  as  a  function  of  0  for  at  an  altitude  of  200  km.  Tbe 
figure  includes  calculated  values  for  o=sig,  o’=sigp,  and  ce=alp.  Then,  the  Schaff  and  Chambre  formulas 
are  evaluated  for  Q  and  Q.  Finally,  tiie  Schamberg  formula  for  CJ  is  shown.  This  last  formula  does  not 
include  any  effects  of  the  Maxwell  thermal  motion,  though  it  uses  tiie  HSN  derived  value  for  a  consistent 
with  c  and  o'.  There  is  a  large  difference  between  the  two  models  for  Q.  Though  it  is  not  shown,  the 
difference  between  the  two  models  for  Q  is  similar. 


Figure  12.  Flat  plate  C^for  Nj. 


A  second  example  is  given  for  He  at  800-km  altitude  in  Figure  13.  Again,  tbe  calculated 
accommodation  (X)efficients  and  the  resulting  values  for  and  C,  are  seen.  In  both  cases,  the  dificrence 
between  the  Schaff  and  Chambre  formulation  and  die  Schamberg  formulation  is  substantial,  often 
exceeding  a  factor  of  two.  In  both  cases,  the  Schaaf  and  Chambre  fonnulation  gives  different  results  fiom 
the  Sdtamberg  approach. 


83  C4  FOR  SPHERES 

The  calculation  of  for  spherical  satellites  can  now  be  compared.  In  this  case,  C,,  is  calculated 
for  each  thermosphere  constituent  at  a  nuntber  of  altitudes.  It  is  calculated  using  die  Schaaf  and  Chambre 
model  the  Schambt^g  model  as  modified  by  Alfonso  et  al.  Cl,  and  the  Schamberg,  Alfonso  ct  al. 
model  u^ng  Uk  Goodman  and  Wachtuan  formula  for  a  CJ\  CJ  attempts  to  account  for  the  thermal 


motion  of  tile  Maxwell  gas  and  the  variation  of  a  witii  0,  and  Q"  only  accounts  for  the  thermal  motion. 
Recall  that  the  quadrature  in  6  runs  from  -id!  to  %I2.  The  HSN  model  evaluated  for  6<0  gives  c=c'~aM). 


Figure  13.  Flute  plate  Cifor  helium. 


There  are  important  contributions  to  for  6<0,  especially  when  5=1,  This  is  the  case  for  the  last 
exan^le  given;  there,  the  thermal  tenqierature  is  assunied  to  be  10000  K,  This  last  case  is  the  one  studied 
by  Alfonso  et  al.,  the  Lageos  satellite,  that  exhibits  anomalous  acceleration. 

Table  13  gives  the  detailed  results.  In  Figure  14,  Q  is  plotted  as  a  function  of  velocity  for  some 
of  the  thermosphere  constituents.  In  general  the  agreement  between  the  Scbaaf  and  Chambre  model  and 
the  Schamberg,  Alfonso  et  al.  model  is  always  less  than  ±10%.  However,  the  disagreement  between  the 
Scbaaf  and  Chambre  model  and  the  Schamberg  model  with  constant  a  can  be  as  much  as  21%,  and  is 
generally  greater  than  10%.  This  is  surprising  given  the  much  larger  disagreement  when  the  flat  plate 
models  are  compared. 
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Figure  14,  Q  for  spherical  satellites. 


Now  for  altitades  less  than  1000  km,  the  agreement  of  the  three  models  for  N  and  0  is  2  to  3%. 
However,  for  all  ottier  constituents,  the  difference  ranges  from  1 1  to  19%.  In  particular,  at  800  km,  where 
He  and  0  are  tile  principal  constituents,  the  errors  are  17  and  4%,  respectively. 

For  the  case  of  a  deep-space  sphere,  with  high  Maxwell  ten^atuie,  the  principal  constituent  is 
H.  In  this  case,  the  error  in  C/  is  ■»-21%,  le.,  one  overestimates  the  drag  force  by  21%. 
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TABLE  13 

C(j  for  Spherical  Satellites 


q 

V 

T 

Tw 

M 

Y 

Cd 

Cd' 

Cd" 

e' 

8" 

km 

km/s 

K 

K 

amu 

H 

200 

7.784 

700 

295 

1.008 

1.667 

2.7298 

2.8629 

3.1430 

5 

15 

He 

200 

7.784 

700 

295 

4.003 

1.667 

2.3960 

2.5645 

2.6608 

7 

11 

N 

200 

7.784 

700 

295 

14.004 

1.667 

2.2589 

2.4033 

2.3075 

6 

2 

0 

200 

7.784 

700 

295 

15.999 

1.667 

2.2403 

2.3765 

2.2993 

6 

3 

N2 

200 

7.784 

700 

295 

28.013 

1.400 

2.1265 

2.2049 

2.3741 

4 

12 

02 

200 

7.784 

700 

295 

31.999 

1.400 

2.0913 

2.1447 

2.407 

3 

15 

Ar 

200 

7.784 

700 

295 

39.948 

1.400 

2.0773 

2.0708 

2.4654 

■0 

19 

H 

400 

7.669 

700 

295 

1.008 

1.667 

2.7271 

2.8655 

3.1534 

5 

16 

He 

400 

7.669 

700 

295 

4.003 

1.667 

2.3840 

2.5610 

2.6632 

7 

12 

N 

400 

7.669 

700 

295 

14.004 

1.667 

2.2459 

2.3970 

2.3082 

7 

3 

0 

400 

7.669 

700 

295 

15.999 

1.667 

2.2277 

2.3698 

2.2999 

6 

3 

N2 

400 

7.669 

700 

295 

28.013 

1.400 

2.1180 

2.1963 

2.3744 

4 

12 

O2 

400 

7.669 

700 

295 

31.999 

1.400 

2.0859 

2.1362 

2.4074 

2 

15 

Ar 

400 

7.669 

700 

295 

39.948 

1.400 

2.0796 

2.0705 

2.465 

-0 

19 

H 

600 

7.558 

700 

295 

1.008 

1.667 

2.7252 

2.8681 

3.1637 

5 

16 

He 

600 

7.558 

700 

295 

4.003 

1.667 

2.3727 

2.5576 

2.6657 

8 

12 

N 

600 

7.558 

700 

295 

14.004 

1.667 

2.2337 

2.3908 

2.3088 

7 

3 

0 

600 

7.558 

700 

295 

15.999 

1.667 

2.2159 

2.3632 

2.3004 

7 

4 

N2 

600 

7.558 

700 

295 

28.013 

1.400 

2.1105 

2.1879 

2.3748 

4 

13 

0 

600 

7.558 

700 

295 

31.999 

1.400 

2.0817 

2.1282 

2.4077 

2 

16 

At 

600 

7.558 

700 

295 

39.948 

1.400 

2.0810 

2.0704 

2.4658 

-1 

18 

H 

800 

7.452 

700 

295 

1.008 

1.667 

2.7240 

2.8707 

3.1741 

5 

17 

He 

800 

7.452 

700 

295 

4.003 

1.667 

2.3621 

2.5543 

2.6682 

8 

13 

N 

800 

7.452 

700 

295 

14.004 

1.667 

2.2222 

2.3847 

2.3094 

7 

4 

0 

800 

7.452 

700 

295 

15.999 

1.667 

2.2048 

2.3568 

2.3010 

7 

4 

N2 

800 

7.452 

700 

295 

28.013 

1.400 

2.1039 

2.1797 

2.3751 

4 

13 

O2 

800 

7.452 

700 

295 

31.999 

1.400 

2.0788 

2.1205 

2.4080 

2 

16 

Ar 

800 

7.452 

700 

295 

39.948 

1.400 

2.0819 

2.0702 

2.4661 

-1 

18 

H 

1000 

7.350 

700 

295 

1.008 

1.667 

2.7820 

2.8738 

3.1844 

3 

14 

He 

1000 

7.350 

700 

295 

4.003 

1.667 

2.4102 

2.5511 

2.6707 

6 

11 

N 

1000 

7.350 

700 

295 

14.004 

1.667 

2.2574 

2.3787 

2.3101 

5 

2 

0 

1000 

7.350 

700 

295 

15.999 

1.667 

2.2371 

2.3504 

2.3015 

5 

3 

N2 

1000 

7.350 

700 

295 

28.013 

1.400 

2.1164 

2.1717 

2.3754 

3 

12 

O2 

1000 

7.350 

700 

295 

31.999 

1.400 

2.0858 

2.1134 

2.4083 

1 

15 

Ar 

1000 

7.350 

700 

295 

39.948 

1.400 

2.0826 

2.0700 

2.4663 

-1 

18 

H 

5800 

5.721 

10000 

295 

1.008 

1.667 

6.1588 

5,8562 

7.4522 

-5 

21 

He 

5800 

5.721 

10000 

295 

4.003 

1.667 

4.2973 

3.9866 

4.4837 

.7 

4 

N 

5800 

5.721 

10000 

295 

14.004 

1.667 

3.1747 

2.8608 

2.8196 

-10 

-11 

0 

5800 

5.721 

10000 

295 

15.999 

1.667 

3.0663 

2.7774 

2.7487 

.9 

-10 

N2 

5800 

5.721 

10000 

295 

28.013 

1.400 

2.5311 

2.3946 

2.6439 

-5 

4 

O2 

5800 

5.721 

10000 

295 

31.999 

1.400 

2.3808 

2.2974 

2.6473 

4 

11 

Ar 

5800 

5.721 

10000 

295 

39.948 

1.400 

2.2660 

2.2118 

2.6632 

-2 

18 

Causing  Schaaf  and  Cliainbie  fonnulas  \rith  variable  o(6)  and  a'(0). 

CJ  using  Schamberg,  Alfonso  et  al.  with  variable  a(8). 

C/ using  Schamberg,  Alfonso  et  al.  with  constant  a  fiom  Goodman  and  Wachmaa 

6‘-100(Crf'-C^/Q. 

8“-lOO(C/-Crfi/Crf 
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Figtces  IS  tbiough  21  present  the  dq>endence  of  C£f^  for  each  thomiosph^  coustitaent  as  a  function  of 
vdocity  for  a  number  of  temperatures.  Figures  22  through  2S  presmit  the  dqiendence  of  for  each  thomosphere 
constituent  as  a  function  of  ten^etature  for  a  number  of  velocities. 


8.4  COMPARISON  OF  C*  WITH  HERRERO 

Hararo  [30,31]  tried  to  calculate  tbe  effective  for  a  cylindrical  spacea  .ft  with  motion  along 
the  axis  of  symmetry.  He  showed  tiiat  die  contribution  to  Q  of  the  lateral  surface,  parallel  to  the  motion, 
is  2(Ur)Cis,  where  L  is  the  length  of  the  cylinder,  r  is  the  radius  of  die  cylinde.t,  and  is  the  drag 
coefficient  for  the  surface  element.  Thou^  a  novel  jqjproach,  he  estimated  that  0.07<Cij<0.06.  The 
theory  here  gives,  for  at  200-km  altitude,  Cu=0.0378.  If  we  adopt  Hcirero's  estimate  for  the 
contribution  of  tbe  nose  cone  as  C^pTS,  and  I/r=10,  we  obtain 


C^=Ccy+2yCj^=lJ+0.76=2.26 
Herrero's  estimate  ranged  between  2.7  and  2.9. 


23174M5 


VBIOCITY  (ton/s) 


Figure  15.  Cj  for  sphere:  hydrogen. 
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2S1741-1S 


6  7  8 

VELOCITY  (km/s) 

Figure  16.  Q  for  sphere:  helium. 


VELOCITY  (km/t) 


Figure  17,  for  sphere:  nitrogen. 
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VELOCITY  (km/s) 

Figure  20.  Q  for  sphere:  0^. 


231741-21 


VELOCITY  (kiiV») 

Figure  21.  Q  for  sphere:  argoru 
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231741.22 


I'EMPERATURE  (K) 


Figure  24.  Q  for  sphere:  nitrogen. 


Figure  25.  for  sphere:  oxygen. 


^9 


9.  SUMMARY 


We  have  developed  a  model  for  scattering  that  is  appropriate  for  satellite  conditions,  i.e.,  satellite 
velocities  and  materials  and  thermosphere  constituents  and  temperatures.  The  model  is  based  on  the 
formalism  of  Hurlbut,  Sherman,  and  Nocilla  (the  HSN  model)  and  uses  published  data  on  momentum 
accommodation  coefficients  to  determine  the  underlying  parameters.  The  model  agrees  with  observed 
momentum  accommodation  coefficients  to  2.5%.  Though  applied  to  the  seven  atmospheric  constituents, 
the  model  is  limited  to  velocities  between  1.77  and  13.66  km/sec  and  to  scattering  surfaces  that  can 
reasonably  be  viewed  as  presenting  oxygen  atoms  to  the  impinging  particles. 

The  HSN  model  predicts  the  three  accommodation  coefficients,  a,  -5',  and  a,  as  a  function  of 
incident  angle.  This  allows  a  comparison  of  the  ballistic  coefficient  model  of  Schaaf  and  Chambre  [15] 
(based  on  a  and  o’)  and  that  of  Schamberg  [10,11]  (based  on  a).  Schamberg's  model  is  the  basis  of 
many  aeronomy  analyses  [4].  Schamberg's  model  does  not  account  for  the  thermal  motion  of  the 
impinging  molecules  and  assumes  the  scattered  molecules  are  fully  accommodated  and  scattered  diffusely. 
Alfonso  et  al.  [16]  tried  to  extend  Schamberg’s  model  to  account  for  atmospheric  thermal  motion  for 
spherical  satellites.  We  have  been  unable  to  fuUy  reconcile  these  two  {^roaches.  The  Schaaf  and 
Chambre  approach  is  adopted  as  more  fundamental,  making  fewer  assumptions. 

Ihe  agreement  of  the  Schaaf  and  Chambie  model  and  the  Schamberg  model,  using  the  HSN 
theory  to  determine  scattering  properties,  is  no  better  than  a  f  xtor  of  two.  ^en  applied  to  spheres, 
taking  into  account  tiie  aspect  dependence  of  a,  the  two  models  agree  to  better  than  10%.  However,  using 
a  constant  a,  the  Schamberg,  Alfonso  ot  al.  model  has  errors  exceeding  20%. 

For  application  of  these  models  to  calculation  of  satellite  drag,  one  can  expect  an  error  exc^ding 
20%,  using  the  simple  ballistic  coefficient.  Ihe  errors  are  not  constant,  depending  on  altitude  and 
thermosphere  constituent.  Therefore,  one  can  expect  a  significant  contribution  to  the  calculation  of 
satellite  drag  from  unmodeled  errors  in  accommodation  coefficients  and  ballistic  coefficients,  even  for 
spherical  satellites. 
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